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ABSTRACT 

We present detailed predictions for the confusion noise due to extragalactic sources 
in the far-IR/(sub)-millimeter channels of ESA/ISO, NASA/Spitzer, ESA/Herschel 
^vq \ and ES A/Planck satellites, including the contribution from clustering of unresolved 

J> ■ SCUBA galaxies. Clustering is found to increase the confusion noise, compared to the 

Q\ ' case of purely Poisson fluctuations, by 10-15% for the lowest frequency (i.e. lowest an- 

Q\ , gular resolution) Spitzer and Herschel channels, by 25-35% for the 175 /Ltm ISOPHOT 

channel, and to dominate in the case of Planck/HFI channels at v > 143 GHz. Al- 
though our calculations make use of a specific evolutionary model (Granato et al. 
2004), the results are strongly constrained by the observed counts and by data on the 
redshift distribution of SCUBA sources, and therefore are not expected to be heavily 
model dependent. The main uncertainty arises from the poor observational definition 
. of the source clustering properties. Two models have been used for the latter: a power- 

On' law with constant slope and a redshift-independent comoving correlation length, ro, 

and the standard theoretical model for clustering evolution in a ACDM universe, with 
^ ■ a redshift-dependent bias factor. In both cases, the clustering amplitude has been nor- 

^ ' malized to yield a unit angular correlation function at 9q = l"-2" for 850/j,m sources 

■ fainter than 2 mJy, consistent with the results by Peacock et al. (2000). This nor- 

malization yields, for the first model, ro ~ 8.3 Mpc/h, and, for the second model, an 
effective mass of dark matter haloes in which these sources reside of Mhaio ~ 1.8 x 10 13 
M Q /h. These results are consistent with independent estimates for SCUBA galaxies 
',— ] ■ and for other, likely related, sources. 
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1 INTRODUCTION 

Our knowledge of the early galaxy evolution has significantly 
improved in the last years thanks to the observations per- 
formed in the sub-millimeter waveband, in particular with 
the 450-850/im Sub-millimeter Common User Bolometric 
Array (SCUBA) camera (Holland et al. 1999). The sources 
detected by SCUBA likely represent the progenitors of the 
local giant elliptical galaxies whose properties indicate that 
the bulk of their star formation activity occurred at high 
redshift and in a relatively short time (Thomas et al. 2002; 
Romano et al. 2002). SCUBA galaxies are in fact located at 
high (> 2) redshifts (Dunlop 2001; Ivison et al. 2002; Chap- 
man et al. 2002, 2003) and display an enormous star forma- 
tion rate (~1000 Mo/yr), which allows masses in stars of 
the order of 10 11 Mq to be assembled in times shorter than 



1 Gyr (Smail et al. 2002; Aretxaga et al. 2003; Chapman et 
al. 2003; Saracco et al. 2003). 

For dark matter to stellar mass ratios typical of massive 
ellipticals (see, e.g., Marinoni & Hudson 2002; McKay et al. 
2002), the dark matter haloes in which the sub- millimeter 
galaxies reside have masses of > 10 13 M Q . Since these mas- 
sive haloes sample the rare high-density peaks of the pri- 
mordial dark matter distribution (Kaiser 1984; Mo & White 
1996), SCUBA sources are expected to exhibit a strong spa- 
tial clustering, similar to that measured for Extremely Red 
Objects (EROs; Daddi et al. 2001, 2003), for which similar 
masses have been inferred (Moustakas& Somerville 2002). 

Direct measurements of clustering properties of SCUBA 
galaxies are made difficult by the poor statistics and by the 
fact that they are spread over a wide redshift range, so that 
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their clustering signal is strongly diluted. However, tentative 
evidences of strong clustering with a correlation length of ~ 
8-13 Mpc/h, consistent with that found for EROs, have been 
reported (Webb et al. 2003; Smail et al. 2003). Peacock et 
al. (2000, henceforth POO), after having removed all sources 
brighter than 2 mjy, found, in the ~ 2 x 2 arcmin 2 850/im 
map by Hughes et al. (1998) in the Hubble Deep Field, some 
evidence for clustering of the background source population 
consistent with an angular correlation function of the form 
w(6) = (e/Ov)- -* with 0o = l"-2". 

As first pointed out by Scott & White (1999), the clus- 
tering signal may provide the dominant contribution to the 
power spectrum of small-scale fluctuations due to extra- 
galactic sources at (sub)-millimeter wavelengths. Detailed 
calculations were carried out by Haiman & Knox (1999), 
Knox et al. (2001), Magliocchetti et al. (2001), and Perrotta 
et al. (2003). The general conclusion is that the amplitude 
of the clustering signal on sub-degree scales is comparable 
to that from the cosmic microwave background at 850 /im, 
and quickly overwhelms it at shorter wavelengths. 

We address this issue from a different point of view, 
i.e. we analyze the effect of clustering in increasing the 
confusion noise (and therefore the detection limit) in far- 
IR/sub- millimeter surveys. The effect is particularly rele- 
vant for confusion limited surveys with poor spatial reso- 
lution, since, for small enough angular scales, the ratio of 
Poisson-to-clustering fluctuations decreases with increasing 
angular scale (De Zotti et al. 1996). 

In particular, we consider the surveys in the ISOPHOT 
175^m channel, in the Spitzer (formerly SIRTF) MIPS 
160 /jm channel, and those performed with Herschel/SPIRE, 
and with Planck/HFI. Essentially all current estimates of 
the confusion noise for these instruments (Blain et al. 1998; 
Dole et al. 2001; Xu et al. 2001; Franceschini et al. 2001; 
Rowan-Robinson 2001; Lonsdale et al. 2003; Lagache et al. 
2003; Dole et al. 2003; Rodighiero et al. 2003) do in fact only 
take into account the Poisson contribution. An attempt to 
allow also for the effect of clustering was worked out by Tof- 
folatti et al. (1998), with reference to the Planck mission. 

For our calculations we adopt the model by Granato et 
al. (2004, henceforth GDS04), which successfully reproduces 
the observed counts and the available information on the 
redshift distribution of SCUBA sources. The data actually 
tightly constrain the volume emissivity of these sources as a 
function of cosmic time, one of the key ingredients to derive 
the amplitude of the fluctuations. As mentioned above, the 
other key ingredient, namely the correlation function, is only 
poorly known, although some observational hints exist. 

The layout of the paper is as follows. In Section 2 we 
describe the formalism. Section 3 summarizes the main as- 
pects of the model by Granato et al. (2004). In Section 4 
we illustrate the models for the two-point spatial correla- 
tion function of SCUBA sources. In Section 5 we present 
and discuss our main results. 

Throughout the paper we will assume a flat cold 
dark matter (CDM) cosmology with Qa = 0.7 and h = 
.Ho/100 kms -1 Mpc -1 = 0.7, in agreement with the first- 
year WMAP results (Spergel et al. 2003). 



2 CONFUSION NOISE DUE TO CLUSTERED 
SOURCES: THE FORMALISM 

The autocorrelation function, C(Sd,d*), of intensity fluctu- 
ations due to sources fainter than the detection limit Sd, as 
a function of the angular separation , writes (see De Zotti 
et al. 1996): 



C{S d ,6 t ) = j dz- L dz 2 J dUdL^ 

I 

Jb 



K(zi) L 2 K(z 2 ) 



4nd 2 L (zi) 47Td|(2 2 ) 

dn 1 dn 2 {6M(L 1 ,r 1 )6M{L 2 ,r 2 ))f(e 1 ,4> 1 )f{e2,4>2) ■ (1) 



Here dh{z) is the luminosity distance, K(z) = (1 + z)L[v(l + 
z)]/L[is] is the K-correction for monochromatic observations 
at the frequency v, the vectors ri and r 2 represent the 
two directions of observation separated by an angle , and 
8M{L, r) is the fluctuation around the mean number density 
M{L, z) of sources with luminosity L and redshift z: 



Af(L,z) = 



dV 



(z)*(L,z) , 



(2) 



dzdQ. 

dV/dzdQ, being the comoving volume element, and $ the 
luminosity function. The angles (9i, 0i) and (9 2 , 4> 2 ) are 
measured from the axes ri and r 2 , respectively. 

The integration in redshift is made within the in- 
terval Z =[2 m i n , Zmax], where the sources are expected 
to exist. The integration over L runs over the interval 
£ =[L m in, min[L max , L(Sd, z)]], where L max and L min are the 
maximum and minimum intrinsic luminosities, and L(Sd,z) 
corresponds to the flux density Sd at the redshift z. 

The spatial response function of the detector, f(9,(f>), 
is assumed to be axially symmetric and Gaussian: 



/(M)=e- (9/8)/2 , 

where O relates to the beam FWHM through 
FWHM 



2\/2ln2 



(3) 



(4) 



If the luminosities of the sources are statistically indepen- 
dent of their positions, we get (Dautcourt 1977): 

<5M{L 1 ,r 1 )8M{L 2 ,v 2 )> = 

U(L 1 ,z 1 )S D (r 2 -r 1 )S D (L 2 -L 1 ) 
+ AA(ii,zi)AA(Li,z 2 )^(r 2 -r 1 ,z) , (5) 

S D being the Dirac function. The first term on the right- 
hand side corresponds to a Poisson distribution of point 
sources, while the second one is due to clustering with a 
two-point spatial correlation function £. 

The variance of intensity fluctuations within the tele- 
scope beam, whose square root, on, will be referred to as 
confusion noise, is obtained setting = 0: 

o N (S d ) = lim C(S d , 9*) = op + oc , (6) 

B * — »0 

and is the quadratic sum of the Poisson term: 
/ dn/ 2 (0,0)x 

J beam 



2 

CTp 



dL 



LK{z) 

^dW) 



*{L,z), (7) 



and of the clustering contribution, which, in the small angle 
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approximation (i.e. 9\, 62 « 1) and assuming that the max- 
imum scale of appreciable clustering is much smaller than 
the Hubble radius, c/Ho, writes: 

<rc = Cj-jf) 2 / d£lid£km,4>i)f(02,<h) x 

J beam 

J z {i + zy&wJ^ c/h{ z) ^i_ (d#w/r) 2 ' 1 

where r is a comoving scale and j e s(z, So) is the effective 
comoving volume emissivity at redshift z contributed by 
sources with fluxes S < Sd'- 

jes(z,S d ) = J dLK(z)L$(L,z) . (9) 

In Eq. (8) E(z) describes the time evolution of the Hubble 
parameter [H(z) — HqE{z)\, d$(z) is the comoving linear 
distance corresponding, at a given redshift z, to an angular 
separation ■& in the sky*. The upper limit of integration over 
r, r sup , is the comoving scale at which the correlation van- 
ishes; we adopt r sup = 30Mpc/h, i.e. approximately three 
times the measured clustering radius ro (see Section 4). 

The detection limit, Sd, for a confusion limited survey, 
is defined by: 

S d = qx a N (S d ) , (10) 
where the parameter q is usually chosen in the range [3,5]. 

3 THE EVOLUTIONARY MODEL 

According to the GDS04 model, SCUBA sources are inter- 
preted as spheroids observed during their major episode of 
star-formation, whose evolution and duration, shorter for 
more massive objects, is substantially affected not only by 
supernova feedback, but also by the growth of a central 
super- massive black hole (SMBH), and by the ensuing QSO 
activity. The relative importance of the QSO feedback, com- 
pared to supernovae, increases with the galaxy mass. This 
scenario, already explored in a phenomenological way by 
Granato et al. (2001), has been substantiated by GDS04 in 
a model which follows, using simple semi-analytic prescrip- 
tions, the evolution of baryons within dark matter haloes 
with total mass M v i r >2. 5 x 10 11 M Q , forming at z > 1.5 at the 
rate predicted by the canonical hierarchical clustering sce- 
nario. The mass and redshift cuts are meant to crudely filter 
out the haloes associated with disk and irregular galaxies, 
whose formation is not quantitatively addressed by GDS04. 
However, disk and irregular galaxies are envisaged as asso- 
ciated to haloes virializing at lower redshifts, eventually in- 
corporating most of the less massive haloes which virialized 
at earlier times and may become galactic bulges. 

We note that according to GDS04 (see also Granato 
et al. 2001), the high-redshift QSO activity marks and con- 
curs to the end of the major episode of star formation in 
spheroids. Thus, there is a clear evolutionary link between 
SCUBA sources and high-redshift QSOs. Also, in the mass 



* In the small angle approximation: $ 2 = d\ + 9% — 19\Qi cos(0i — 
<t>2) 



and redshift range considered, the model implies a one-to- 
one correspondence between haloes and spheroidal galax- 
ies, consistent with the available data on the clustering of 
Lyman-break galaxies (Bullock et al. 2002) and SCUBA 
galaxies (Magliocchetti et al. 2001) 

The GDS04 model predicts a local SMBH mass function 
and a relationship between black hole and spheroid mass (in 
stars) which is in excellent agreement with the observational 
data. Also, when coupled with the most updated version of 
GRASIL (whose original version is described by Silva et al. 
1998), the code which computes in a self-consistent way the 
chemical and spectro-photometric evolution of galaxies from 
far-UV to radio, it yields predictions which are fully consis- 
tent with a number of observables extremely challenging for 
all other semi-analytic models, including the sub-millimeter 
number counts and redshift distribution of spheroidal galax- 
ies, as well as the epoch-dependent K-band luminosity func- 
tion. 

High-redshift star-forming spheroids are not the only 
important class of objects in the far-infrared and sub- 
millimeter wavebands. We also need to take into account 
late-type, normal galaxies and starburst galaxies which in 
fact dominate the bright tail of the counts. To deal with 
these sources, we adopt, as usual, a phenomenological ap- 
proach (e.g. Franceschini et al. 2001; Takeuchi et al. 2001; 
Rowan-Robinson 2001; Gruppioni et al. 2002; Lagache et 
al. 2003), which consists in simple analytic recipes to evolve 
their 60/im local luminosity functions (Saunders et al. 1990) 
and appropriate templates for the spectral energy distribu- 
tions (SED) to extrapolate the models to different wave- 
lengths. The prescriptions adopted here are those defined 
by Silva et al. (2004), which provide reasonable fits to the 
results of IRAS and ISO surveys. Briefly, normal late-type 
galaxies are assumed to undergo a pure luminosity evolution 
described by L(z) = L(0)(1 + z) 1 ' 5 and to have the SED of 
the Sc galaxy NGC 6946. The SED of the galaxy NGC 6090 
is adopted for starburst galaxies and the luminosity function 
of these objects, $sb, is assumed to evolve both in luminos- 
ity and in density: & SB [L(z),z] = $sb[L(z)/(1 + z) 2 ' 5 , z = 
0](l + z) 3,5 . The evolutionary laws for both populations were 
assumed to apply up to z = 1; then the luminosity functions 
were kept constant up to z cuto B = 1.5. 

Spiral and starburst galaxies are relatively weakly clus- 
tered (Fisher et al. 1994; Loveday et al. 1995; Guzzo et al. 
1997; Madgwick et al. 2003), so that their contributions to 
the confusion noise turn out to be dominated by Poisson 
fluctuations. In the following we will therefore neglect their 
contribution to ac- 

At millimeter wavelengths the composition of the bright 
counts drastically changes since flat spectrum radio sources 
take over. Their Poisson fluctuations have been estimated 
by extrapolating the 15 GHz differential counts, n(Sj y ) = 
5157/ 15 Jy _1 sr _1 (Waldram et al. 2003), with a spectral 
index equal to zero (i.e. assuming a flux density constant 
with frequency), consistent with the WMAP results (Ben- 
nett et al. 2003). This choice for the spectral index is likely 
to increasingly overestimate the counts with increasing fre- 
quency, as the radio-source spectra steepen both because 
emitting regions become optically thin and because of elec- 
tron ageing effects. However, this has little impact on our 
results since at sub-millimeter wavelengths radio sources are 
anyway a minor component. 
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Although radio sources are rather strongly clustered 
(Magliocchetti et al. 1998; Blake & Wall 2002a,b; Overzier et 
al. 2003, Magliocchetti et al. 2004), their clustering signal is 
highly diluted due to the broadness of their luminosity func- 
tion (e.g. Dunlop & Peacock 1990) to the effect that their 
ac is small when compared to ap, and can be neglected 
(Argiieso et al. 2003; Toffolatti et al. 1998). 



4 THE TWO-POINT CORRELATION 
FUNCTION OF SCUBA GALAXIES 

POO carried out a power-spectrum analysis of the 850/im 
map of the northern Hubble Deep Field by Hughes et al. 
(1998), after subtracting sources brighter than Sd = 2mJy. 
They found, for multipole numbers I in the range 10 4 $ 
£ £ 5 x 10 4 , some evidence of power in excess of the sum of 
instrumental noise and of estimated Poisson fluctuations due 
to unresolved sources. Such excess power can be accounted 
for by source clustering described by an angular correlation 
function of the form: 



w{B) = (0/0q) 



(11) 



with 8o in the range 1-2 arcsec. As made clear by POO, 
however, the estimated amplitude of Poisson fluctuations 
is rather uncertain due to our poor knowledge of the 850/im 
counts below 2 mjy, and the possibility that they account 
entirely for the detected confusion noise cannot be ruled out. 
The more recent data on the faint 850 /im counts (Chapman 
et al. 2002; Smail et al. 2002; Knudsen & van der Werf 2003) 
are in good agreement with the analysis by POO and do not 
allow firmer conclusions on the amplitude of the clustering 
signal. On the other hand, as shown below, an angular corre- 
lation function with the amplitude suggested by POO is con- 
sistent with a number of other data on clustering of SCUBA 
galaxies themselves (although, again, the significance of the 
clustering detection is not very high) and of other popula- 
tions which are closely linked to such galaxies, such as EROs. 
We thus conclude that Eq. (11) with do — l"-2" provides 
the current best guesstimate for the normalization of our 
models. 

The spatial correlation function, £ sp h, can be obtained 
by inverting the Limber's (1953) equation: 



w(e,s d ) 



i 



-^850 i^d 



dz 



\4ttH J (1 



dr 



Jeff (z, Sri) 

(l + z^E^z) 
£s P h(r,z) 



dg(z) (c/H(z)) v /l_( deW/r )2 



, (12) 



with Sd ~ 2mJy. The quantity /gso is the 850/im back- 
ground intensity produced by sources fainter than Sd- The 
other symbols have the same meaning as in Eq. (8). As 
mentioned in the Introduction, SCUBA galaxies are inter- 
preted by GDS04 as forming massive spheroidal galaxies. 
The model yields I 85 o(S d ) = 6.6 x 10~ 19 erg/s/cm 2 /Hz/sr. 

In the following, we will consider both a phenomeno- 
logical, power-law model for £ sp h (Model 1) and a physically 
motivated model (Model 2). Both models are normalized to 
Eq. (11) with 6*o = l"-2" in the range of scales probed by 
POO (see Fig. 1). 
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Figure 1. Comparison of the adopted models for the angular 
correlation function with Eq. (11) with 9q = 1"— 2", in the interval 
of angular scales probed by POO. The long-dashed curve has been 
obtained using the prescriptions by Peacock & Dodds (1996) for 
generating the dark matter power spectrum, while the dotted 
curve corresponds to the more accurate method by Smith et al. 
(2003) 



4.1 Model 1 

Our first model assumes the usual power-law shape: 



£s P h(r,z) = [r/r (z)Y 



(13) 



In view of the tight connection between spheroidal galaxies 
and active nuclei at their centers, entailed by the GDS04 
model, we further assume that the correlation length, ro, is 
constant in comoving coordinates, as suggested by quasar 
data (Croom et al. 2001; Outram et al. 2003) at least in the 
redshift interval 0.5 < z < 2.5-3. 

Then, from Eq. (12), we find that the range of values of 
ro corresponding to 9o = l"-2" is: 



r = 8.3 ± 1.3 Mpc/h , 



(14) 



very close to the range of values found for bright 2QZ 
quasars (r = 8.37 ± 1.17Mpc/h; Croom et al. 2002), and 
fully consistent with the tentative estimate of the correlation 
length of SCUBA galaxies obtained by Smail et al. (2003): 
ro = 8 ± 3 Mpc/h. The above value of the comoving corre- 
lation length is also compatible with the ones measured for 
bright local elliptical galaxies (Loveday et al. 1995; Guzzo 
et al. 1997; Norberg et al. 2002; Madgwick et al. 2003), sup- 
porting the link between these objects and the population 
of SCUBA sources. 



4.2 Model 2 

The spatial correlation function of visible galaxies can be 
described as the product of the correlation function of dark 
matter, £DM(r, z), with the bias parameter which relates the 
distribution of light to that of the matter (Matarrese et al. 
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Table 1. Estimated confusion noise and detection limits, Sd (neglecting instrumental noise as well as emissions within our Galaxy and 
fluctuations of the cosmic microwave background) for several far-IR/(sub)-millimcter instruments. Shown arc the Poisson contributions 
(up) due to spiral galaxies (sp), star-forming galaxies (sb), radio-galaxies (rg) and star-forming spheroids (sph). For the latter, we also 
give the clustering fluctuations (<Jc,sph)- At each wavelength, three cases are shown: Poisson contributions only (first line), Poisson plus 
clustering contribution from star-forming spheroids with a correlation function of the form ^ sp h('", z) = (r/ro)~ 1-8 with tq = 8.3Mpc/h 
(Model 1 , second line) , Poisson plus clustering contribution from star-forming spheroids with a correlation function of the form £ sp h (r, z) = 
b 2 (M eB ,z)^- DM with M off = 1.8 X 10 13 M Q /h (Model 2, third line). We have assumed q = 5, so that S d = 5 X (crp sp + up ab + <r p rg + 
CT 2 , 2 u/2 

U P,sph + a C,sph> 
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Figure 2. Confusion noise crjv as a function of the detection limit for spiral galaxies (short-dashed line), starburst galaxies (long- 
dashed line), radio-sources (dotted— long-dashed line) and star-forming spheroids (dotted— short-dashed line), all assumed to be randomly 
distributed. The confusion noise due to clustered star-forming spheroids (Model 1 with ro = 8.3 Mpc/h) is shown by the solid line. The 
vertical dotted lines correspond to the detection limit Sd obtained by taking into account both Poisson and clustering fluctuations. 



1997; Moscardini et al. 1998). £dm(t, z) is determined by 
the cosmology and by the primordial density perturbations 
(Peacock & Dodds 1996; Smith et al. 2003), for which we 
have adopted a CDM spectrum, with an index n = 1, a 
shape parameter V = 0.2 and a normalization as = 0.8 
(see e.g. Lahav et al. 2002, Spergel et al. 2003). We will use 
the practical fitting formula by Peacock & Dodds (1996) 
to evolve the dark matter power spectrum into the non- 
linear regime. For the cosmological model adopted here, the 
angular correlation function obtained using the improved 
treatment by Smith et al (2003) is only slightly different 
and is intermediate between those yielded by our two models 
(Fig. 1). 

Under the assumption of a single object per halo, which 
is expected to hold for massive objects at high redshift (z £ 
1) and is implicit in the GDS04 model, the bias parameter, 
b, can be written as a function of an effective dark-matter 
halo mass, M B g (which is close to the minimum halo mass, 
cf. Moscardini et al. 1998), and of redshift (see, e.g., Mo & 



White 1996; Sheth & Tormen 1999). We can then write the 
spatial correlation function of star-forming spheroids as 

£ sp h(r, z) = b 2 (M cS , z)£ D M(r, z) . (15) 

The effective mass M ett - is determined by the condition that, 
in the range of angular scales probed by POO, the model 
correlation function matches the measured one. Using the 
analytic formulae by Sheth & Tormen (1999) for the bias 
factor and by Peacock & Dodds (1996) for £dm(t, z) we find: 

M cB = 1.8 x 10 13 A/ o /h , (16) 

consistent with the estimates by Moustakas & Somerville 
(2002) for ERO halo masses (the minimum mass of ERO 
dark matter haloes is estimated to be ~ lO 13 M0/h, while 
the galaxy number-weighted average halo mass is ~ 5 x 
lO 13 M /h). We recall that EROs are either massive dusty 
galaxies like SCUBA galaxies, or evolved giant ellipticals, a 
later evolutionary phase of SCUBA galaxies, in the scenario 
discussed by GDS04. 

For a virialization redshift z v ir ~ 4, GDS04 predict a 



© 0000 RAS, MNRAS 000, 000-000 



Confusion noise at far-IR to millimeter wavelengths 7 



ratio M c ff/M sp h ~ 40, where M sp h is the present-day mass 
in stars of spheroids (see Figure 5 of their paper). Thus M e ff 
corresponds to a mass in stars of ~ 6 x 10 11 Mg, compatible 
with values for the stellar masses inferred for SCUBA galax- 
ies (Smail et al. 2003) and for high-2 galaxies with very red 
near-IR colours (Saracco et al. 2003). 

As a last remark, we note that (as also illustrated by 
Fig. 1) Models 1 and 2 widely differ on small angular scales. 
The small scale flattening of w(6) implied by Model 2 is due 
to the fact that at high redshift (z>2) density perturbations 
are close to the linear regime, even on very small scales. 



5 RESULTS AND DISCUSSION 

Table 1 gives our estimated contributions to the confusion 
noise and the corrisponding detection limits Sd for all the 
populations previously described (obtained by adopting q = 
5 in Eq. (10)), in the case of: 

i) the 175 /im channel of the Imaging Photo-Polarimeter 
of the ESA's ISO satellite (ISOPHOT); 

ii) the longest wavelength channel (160 (im) of the Multi- 
band Imaging Photometer (MIPS) of NASA's Spitzer (for- 
merly SIRTF) satellite, launched last August; 

iii) all channels (250, 350 and 500 fj,m) of the Spectral and 
Photometric Imaging REceiver (SPIRE) of the ESA's Her- 
schel satellite scheduled for launch in 2007; 

iv) all channels (350, 550, 850, 1380, 2100, and 3000 £tm) of 
the High Frequency Instrument (HFI) of the ESA's Planck 
satellite, to be launched jointly with Herschel. 

The values for the angular resolutions (FWHM) have 
been taken from Dole et al. (2001; ISOPHOT), Lons- 
dale et al. (2003; MlPS/Spitzer), Griffin et al. (2000; 
SPIRE/Herschel), and Lamarre et al. (2003; HFI/Planck). 

As discussed above, forming spheroidal galaxies are the 
only population, among those considered here, whose fluctu- 
ations are dominated by clustering. We have therefore only 
presented in Table 1 results for the instruments and the 
wavelengths where this population gives an important con- 
tribution to the confusion noise. For each channel, the first 
line corresponds to pure Poisson fluctuations while the sec- 
ond and the third lines include the contributions of cluster- 
ing based on Model 1 and 2, respectively. The increment of 
Poisson fluctuations when we allow for the effect of cluster- 
ing is obviously due to the increment of Sd- 

The relative importance of contributions from each class 
of objects as a function of the flux limit is illustrated in 
Fig. 2 by the function on [see Eq. (6)]. The flattening of the 
function oc, sp h{Sd) at bright flux densities follows from the 
fact that the main contribution does not come from sources 
just below the detection limit, as is frequently the case for 
Poisson fluctuations, but from sources at redshifts where the 
effective volume emissivity, j e a [see Eq. (8)], is maximum. 

For power-law differential counts, n(S) oc S~ f3 , op oc 
Sf"^^ if P < 3. This power-law behaviour, obeyed by 
op iTg because of our adoption of the power- law representa- 
tion of the counts by Waldram et al. (2003), must break 
down at faint flux densities, where counts must converge, 
and must flatten at bright flux densities (S k, Uy), where 



the slope of the counts approaches the Euclidean value (Ben- 
nett et al. 2003). Since the survey by Waldram et al. (2003) 
covers the flux density range 10-1000 mjy, we may have 
somewhat overestimated fluctuations due to radio sources, 
if the counts start converging not far below 10 mjy. 

The slope of the bright counts of normal late- type galax- 
ies is close to the Euclidean value (/3 ~ 2.5), while that of 
starburst galaxies is somewhat steeper, due to their rela- 
tively strong evolution. Correspondingly, for both popula- 
tions on is a flatter function of Sd than in the case of radio 
galaxies, and is flatter for starburst than for normal galaxies. 

On the other hand, the bright tail of the counts of form- 
ing spheroidal galaxies is extremely steep (/3 > 3), as a con- 
sequence of the combined effect of strong evolution and nega- 
tive K-correction. Thus, the main contribution to o PySph (Sd) 
comes from relatively faint sources and this quantity is es- 
sentially constant for large enough values of Sd- 

Clustering accounts for 10-20% of the total confusion 
noise, depending on the assumed model for £ sp h, for the 
160 fim MIPS channel and the SPIRE channels, but its 
contribution increases up to 35-40% of the total confusion 
noise for the 175 /im ISOPHOT channel, and dominates 
the fluctuations due to extragalactic sources in the case of 
Planck/HFI, except in the longest wavelength channel where 
Poisson fluctuations due to radio sources take over. How- 
ever, higher resolution surveys can be used to subtract such 
sources down to flux densities well below the estimated Sd, 
thus decreasing their contribution to fluctuations. We note 
that a similar task is much more difficult in the case of form- 
ing spheroidal galaxies since the main contribution to their 
clustering fluctuations comes from very faint flux densities. 

As noted above, the ratio of clustering-to-Poisson fluc- 
tuations increases with decreasing angular resolution. In fact 
[see Eq. (7)], o P oc 9 while, if £(r) = (r/r )" 7 with 7 = 1.8, 
oc oc 1 ' 6 (De Zotti et al. 1996). It may also be noted that, 
as a consequence of the different dependence of w(6) on an- 
gular scale (see Fig. 1), the values of oc implied by Model 1 
(shown in Fig. 2) exceed those given by Model 2 for higher 
resolution surveys, while the opposite happens for lower res- 
olution surveys, such as those of Planck/HFI. 

Our estimates - obtained with the inclusion of the clus- 
tering contribution - of the confusion noise for the 175 /im 
ISOPHOT channel are in good agreement with the observa- 
tional determination by Dole et al. (2001) in the FIRBACK 
survey fields, once we allow for the different flux limit. These 
authors find a confusion noise of 45mJy when adopting a 
3cr at detection limit of 135 mjy, while we adopt a 5on limit 
of 327-350 mjy (see Table 1, bearing in mind that o in- 
creases with Sd)- The possibility of a significant contribu- 
tion from clustering in the 175 /im FIRBACK survey was 
discussed by Perrotta et al. (2003) and Dole et al. (2003). 

In spite of the different models used and of the different 
criteria adopted to define the limiting flux densities, our es- 
timates of Poisson fluctuations are in reasonable agreement 
with those by Rowan- Robinson (2001), Dole et al. (2003) 
and Xu et al. (2003) for MlPS/Spitzer and by Lagache et al. 
(2003) for SPIRE/Herschel and HFI/Planck. At the wave- 
lengths (160 and 850 (im) where the models are tightly con- 
strained by the data, differences are generally within 20-30 
per cent, our predictions being on the high side. At other 
wavelengths, our estimates tend to be higher than those of 
other authors by up to 50 per cent. 
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A preliminary attempt to allow for the effect of cluster- 
ing of evolving dusty galaxies was carried out by Toffolatti 
et al. (1998), who assumed f (r) = [r/ra(z)]~ ls and constant 
clustering in physical coordinates (stable clustering model: 
ro(z) = ro(z = 0)(1 + z) 1 ^ 3 ^ 1 ' 8 , in comoving coordinates). 
They further adopted an Einstein-de Sitter universe and set 
ro(z = 0) = 10 Mpc/h. Their model yields a comoving cor- 
relation scale length at z = 3, r (z = 3) ~ 4 Mpc/h, i.e. 
about a factor of 2 smaller than indicated by current data 
(see Sect. 4). Correspondingly, they have probably substan- 
tially underestimated the clustering fluctuations. 

Recently, Takeuchi & Ishii (2003) have estimated the ef- 
fect of clustering with reference to the forthcoming ASTRO- 
F surveys. They find that, between 60 and 170 /im, cluster- 
ing increases by ~ 10 per cent the detection limit derived 
by only considering Poisson fluctuations. Their result is thus 
compatible with our estimate for the MlPS/Spitzer 160/im 
channel. 

Finally, it must be stressed that we have only consid- 
ered fluctuations due to extragalactic point sources. Addi- 
tional contributions of varying importance (depending on 
the wavelength and on the position in the sky of the surveyed 
area) come from Galactic (synchrotron, free-free, interstellar 
dust) and zodiacal emissions. At millimeter wavelengths the 
graininess of the sky is actually dominated by fluctuations 
of the cosmic microwave background. Further fluctuations 
may be produced by Sunyaev-Zel'dovich effects in groups 
and clusters of galaxies. Last, but not least, to determine 
the sensitivity of a survey we also have to allow for instru- 
mental noise. Therefore, the values of S d quoted in Table 1 
must be regarded as lower limits. 
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